Insertion and deletion mutations preserved in SARS-CoV-2 variants.

The insertion/deletion (indel) mutation profiles of SARS-CoV-2 variants, including Omicron, remain unclear. We compared whole-genome sequences from various lineages and used preserved indels to infer the ancestral relationships between different lineages. Thirteen indel patterns from twelve sites were seen in ≥ 2 sequences; six of these sites were located in the N-terminal domain of the viral spike gene. Preserved indels in the coding regions were also identified in the non-structural protein 3 (Nsp3), Nsp6, and nucleocapsid genes. Seven of the thirteen indel patterns were specific to the Omicron variants, four of which were observed in BA.1, making it the most mutated variant. Other preserved indels observed in the Omicron variants were also seen in Alpha and/or Gamma, but not Delta, suggesting that Omicron is phylogenetically more proximal to Alpha. We demonstrated distinct profiles of preserved indels among SARS-CoV-2 variants and sublineages, suggesting the importance of indels in viral evolution.

Trinucleotide Substitutions at Two Locations in the SARS-CoV-2 Nucleocapsid (N) Gene.

The genomes of sarbecoviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), incorporate mutations with short sequence exchanges based on unknown processes. Currently, the presence of such short-sequence exchanges among the genomes of different SARS-CoV-2 lineages remains uncertain. In the present study, multiple SARS-CoV-2 genome sequences from different clades or sublineages were collected from an international mass sequence database and compared to identify the presence of short sequence exchanges. Initial screening with multiple sequence alignments identified two locations with trinucleotide substitutions, both in the nucleocapsid (N) gene. The first exchange from 5'-GAT-3' to 5'-CTA-3' at nucleotide positions 28,280-28,282 resulted in a change in the amino acid from aspartic acid (D) to leucine (L), which was predominant in clade GRY (Alpha). The second exchange from 5'-GGG-3' to 5'-AAC-3' at nucleotide positions 28,881-28,883 resulted in an amino acid change from arginine and glycine (RG) to lysine and arginine (KR), which was predominant in GR (Gamma), GRY (Alpha), and GRA (Omicron). Both trinucleotide substitutions occurred before June 2020. The sequence identity rate between these lineages suggests that coincidental succession of single-nucleotide substitutions is unlikely. Basic local alignment search tool sequence search revealed the absence of intermediating mutations based on single-base substitutions or overlapping indels before the emergence of these trinucleotide substitutions. These findings suggest that trinucleotide substitutions could have developed via an en bloc exchange. In summary, trinucleotide substitutions at two locations in the SARS-CoV-2 N gene were identified. This mutation may provide insights into the evolution of SARS-CoV-2.

Variable number tandem repeats of a 9-base insertion in the N-terminal domain of severe acute respiratory syndrome coronavirus 2 spike gene.

Introduction: The world is still struggling against the pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), in 2022. The pandemic has been facilitated by the intermittent emergence of variant strains, which has been explained and classified mainly by the patterns of point mutations of the spike (S) gene. However, the profiles of insertions/deletions (indels) in SARS-CoV-2 genomes during the pandemic remain largely unevaluated yet. Methods: In this study, we first screened for the genome regions of polymorphic indel sites by performing multiple sequence alignment; then, NCBI BLAST search and GISAID database search were performed to comprehensively investigate the indel profiles at the polymorphic indel hotspot and elucidate the emergence and spread of the indels in time and geographical distribution. Results: A polymorphic indel hotspot was identified in the N-terminal domain of the S gene at approximately 22,200 nucleotide position, corresponding to 210-215 amino acid positions of SARS-CoV-2 S protein. This polymorphic hotspot was comprised of adjacent 3-base deletion (5'-ATT-3'; Spike_N211del) and 9-base insertion (5'-AGCCAGAAG-3'; Spike_ins214EPE). By performing NCBI BLAST search and GISAID database search, we identified several types of tandem repeats of the 9-base insertion, creating an 18-base insertion (Spike_ins214EPEEPE, Spike_ins214EPDEPE). The results of the searches suggested that the two-cycle tandem repeats of the 9-base insertion were created in November 2021 in Central Europe, whereas the emergence of the original one-cycle 9-base insertion (Spike_ins214EPE) would date back to the middle of 2020 and was away from the Central Europe. The identified 18-base insertions based on 2-cycle tandem repeat of the 9-base insertion were collected between November 2021 and April 2022, suggesting that these mutations could not survive and have been already eliminated. Discussion: The GISAID database search implied that this polymorphic indel hotspot to be with one of the highest tolerability for incorporating indels in SARS-CoV-2 S gene. In summary, the present study identified a variable number of tandem repeat of 9-base insertion in the N-terminal domain of SARS-CoV-2 S gene, and the repeat could have occurred at different time from the insertion of the original 9-base insertion.

Insertion/deletion hotspots in the Nsp2, Nsp3, S1, and ORF8 genes of SARS-related coronaviruses.

The genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains many insertions/deletions (indels) from the genomes of other SARS-related coronaviruses. Some of the identified indels have recently reported to involve relatively long segments of 10-300 consecutive bases and with diverse RNA sequences around gaps between virus species, both of which are different characteristics from the classical shorter in-frame indels. These non-classical complex indels have been identified in non-structural protein 3 (Nsp3), the S1 domain of the spike (S), and open reading frame 8 (ORF8). To determine whether the occurrence of these non-classical indels in specific genomic regions is ubiquitous among broad species of SARS-related coronaviruses in different animal hosts, the present study compared SARS-related coronaviruses from humans (SARS-CoV and SARS-CoV-2), bats (RaTG13 and Rc-o319), and pangolins (GX-P4L), by performing multiple sequence alignment. As a result, indel hotspots with diverse RNA sequences of different lengths between the viruses were confirmed in the Nsp2 gene (approximately 2500-2600 base positions in the overall 29,900 bases), Nsp3 gene (approximately 3000-3300 and 3800-3900 base positions), N-terminal domain of the spike protein (21,500-22,500 base positions), and ORF8 gene (27,800-28,200 base positions). Abnormally high rate of point mutations and complex indels in these regions suggest that the occurrence of mutations in these hotspots may be selectively neutral or even benefit the survival of the viruses. The presence of such indel hotspots has not been reported in different human SARS-CoV-2 strains in the last 2 years, suggesting a lower rate of indels in human SARS-CoV-2. Future studies to elucidate the mechanisms enabling the frequent development of long and complex indels in specific genomic regions of SARS-related coronaviruses would offer deeper insights into the process of viral evolution.